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Abstract 

It is commonly stated that we have entered the era of precision cosmol- 
ogy in which a number of important observations have reached a degree of 
precision, and a level of agreement with theory, that is comparable with 
many Earth-based physics experiments. One of the consequences is the 
need to examine at what point our usual, well-worn assumption of homo- 
geneity associated to the use of perturbation theory begins to compromise 
the accuracy of our models. It is now a widely accepted fact that the effect 
of the inhomogeneities observed in the Universe cannot be ignored when 
one wants to construct an accurate cosmological model. Well-established 
physics can explain several of the observed phenomena without introduc- 
ing highly speculative elements, like dark matter, dark energy, exponential 
expansion at densities never attained in any experiment (i.e. inflation), 
and the like. Two main classes of methods are currently used to deal 
with these issues. Averaging, sometimes linked to fitting procedures a 
la Stoegger and Ellis, provide us with one promising way of solving the 
problem. Another approach is the use of exact inhomogeneous solutions 
of General Relativity. This will be developed here. 

1 Introduction 

The observation of its structures show that our Universe is not homogeneous. 
We see voids, groups of galaxies, clusters, superclusters, walls, filaments, etc. 
However, it is usually argued in the literature that the Universe should be 
nearly homogeneous at large scales which is supposed to validate the use of 
Friedmannian models. But how large these scales are and what nearly does 
imply is never precisely stated. 

It has however become, during the last few years, as a widely accepted fact, 
that the effect of the inhomogeneities cannot be ignored when one wants to 
construct an accurate cosmological model up to the regions where structures 
start forming and their evolution becomes non-linear. Three different methods 
have been proposed to deal with this issue: 

1. Linear perturbation theory. However, this method is only valid when both 
the curvature and density contrasts remain small, which is not the case 
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in the non-linear regime of structure formation and where the SNe la are 
observed. 

2. Averaging methods 'a la Buchert', promising, but needing to be improved 
(see pQ and references therein). 

3. Exact inhomogeneous solutions, valid at all scales and exact perturbations 
of the Friedmann background which they can reproduce as a limit with 
any precision. 

The use of such exact solutions shows that well-established physics can ex- 
plain several of the phenomena observed in astrophysics and cosmology without 
introducing highly speculative elements, like dark matter, dark energy, expo- 
nential expansion at densities never attained in any experiment (inflation), and 
the like. Here, we foccuss on the application of a couple of exact solutions of 
general relativity to structure formation and evolution and to the reproduction 
of cosmological data. 

2 Use of exact inhomogeneous models in astro- 
physics and cosmology 

Very few exact inhomogeneous solutions of Einstein's equations have been used 
for these purposes. Those which appear most often are: 

1. The Lemaitre - Tolman (L-T) models [3] which are spherically sym- 
metric dust solutions of Einstein's equations. They are determined by one 
coordinate choice and two free functions of the radial coordinate r cho- 
sen among three independent ones: the energy per unit mass , E(r), of 
the particles contained within the comoving spherical shell at a given r, 
the gravitational mass, M(r), contained in that shell and the Bang time 
function, <s(r), meaning that the Big Bang occurred at different times at 
different r values. The homogeneous FLRW model is one sub-case. 

2. The Lemaitre model j2] (usually known as Misner-Sharp [4]) is not an 
explicit solution but a metric determined by a set of two differential equa- 
tions. It represents a spherically symmetric perfect fluid with pressure 
gradient. Its solution is obtained by numerical integration. 

3. The Quasi-spherical Szckeres (QSS) models [5] are dust solutions of Ein- 
stein's equations with no symmetry at all. They are defined by one co- 
ordinate choice and five free functions of the radial coordinate. The L-T 
and FLRW models are sub-cases. 

4. The spherically symmetric Stephani models [6] have also been used for 
cosmological purpose. They are exact solutions with homogeneous-energy 
density and inhomogeneous-pressure. 

L-T models have been the most widely used in cosmology since they are the 
most tractable among the few available ones cited above. However, QSS models 
are currently slightly coming into play. 

But caution with L-T models is required since: 

• An origin, or centre of spherical symmetry, occurs at r = r c where R(t, r c ) = 
for all t (here, R is the areal radius). The conditions for a regular centre 
were derived by Mustapha and Hellaby [7j. 



2 



• Shell crossings, where a constant r shell collides with its neighbour, create 
undesirable singularities while the density diverges and changes sign. The 
conditions on the 3 arbitrary functions of the model that ensure none be 
present anywhere in an L-T model are given in [5] . 

• The assumption of central observer, generally retained for simplicity, can 
be considered as grounded on the observed quasi-isotropy of the CMB 
temperature, and thus as a good working approximation at large scales. 
At smaller scales, it gives simplified models of the Universe averaged over 
the angular coordinates around the observer, i. e., with the relax of only 
one degree of symmetry as regards the homogeneity assumption. The 
central observer location is just an artefact of this angular averaging. 

However, models assuming a non-central observer and L-T Swiss-cheeses 
have also been studied to get rid of possible misleading features of spherical 
symmetry. 

3 Structure evolution 

3.1 Structure evolution with L— T and Lemaitre models 

The L-T class of solutions have been used to reproduce the formation and 
evolution of structures from the near homogeneity seen in the CMB temperature. 

In [5] , the evolution from an initial density profile to a galaxy cluster whose 
density profile approximate the 'Universal Profile' of an Abell cluster implies 
a density amplitude at the initial time which differs from the observed values 
at last scattering by three orders of magnitude. However, the velocity ampli- 
tude at the initial time obtained by the same process is on the border of the 
observationally implied range. This shows that velocity perturbations generate 
structures much more efficiently than density perturbations. Moreover, it is 
demonstrated in [7] that a smooth evolution can take an initial condensation 
to a void. This implies that the initial density distribution does not determine 
the final structure which will emerge from it. The velocity distribution can 
obliterate the initial setup. 

However, a void consistent with the observational data (density contrast less 
than 8 — —0.94, smooth edges and high density in the surrounding regions) is 
very hard to obtain with L-T models without shell crossing [TO] , 

But adding a realistic distribution of radiation (using Tolman models) helps 
forming such voids [11] . 

3.2 Structure evolution with quasi-spherical Szekeres mod- 
els 

In the works reported in the above subsection, the spherically symmetric L— 
T and Lemaitre models were used to study structure formation. However, the 
structures observed in the Universe are far from being spherical and the analysis 
needs thus to be refined by considering a wider variety of astrophysical objects 
with different shapes. The investigation of the evolution of small voids inside 
compact clusters and of large voids surrounded by walls or filaments has been 
performed with QSS models. 

A void with an adjourning supercluster evolved inside an homogeneous back- 
ground (FIGURE I) has been considered in [12]. To estimate how two neigh- 
bouring structures influence each other, the evolution of a double structure in 
QSS models (FIGURE 2) has been compared with that of a single structure 
in L-T models and in linear perturbation theory. In the QSS models studied, 
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the growth of the density contrast is 5 times faster than in the corresponding 
L-T models and 8 times faster than in the linear approach. This could imply 
a strong improvement of the phenomenon of structure formation which, in the 
standard CDM model, is too slow to form the structures corresponding to the 
observations made at larger and larger redshifts. However, the evolution of the 
void is slower within the Szekeres model than it is in the L-T model. This 
suggests that single, isolated voids evolve much faster than the ones which are 
in the neighbourhood of large overdensities where the mass of the perturbed 
region is above the background mass. 

The model of triple structure is composed of an overdense region at the 
origin, followed by a small void which spreads to a given r coordinate. At a 
larger distance from the origin, the void is huge and its larger side is adjacent to 
an overdense region [13] (FIGURE 3). Where the void is large, it evolves much 
faster than the underdense region closer to the 'centered' cluster. The exterior 
overdense region close to the void along a large area evolves much faster than 
the more compact supercluster at the centre (FIGURE 4). This confirms that, 
in the Universe, small voids surrounded by large high densities evolve much 
more slowly than large isolated voids. 

4 The dark energy puzzle 

Since its discovery during the late 1990s OUS], the 'dimming' of distant type la 
supernovae has been mostly ascribed to the influence of a mysterious dark energy 
component, i.e. an unknown 'fluid' or 'field' with negative pressure. Formulated 
in a Friedmannian framework, based upon the 'cosmological principle', this 
interpretation has given rise to the 'Concordance' model where the Universe 
expansion is accelerated by the dark energy pressure. 

However, what we observe is not an accelerated expansion (this is only the 
outcome of the Friedmannian assumption) but the dimming of the supernovae 
as regards their luminosity predicted in the previous Einstein-de Sitter standard 
model of cosmology. More exactly we establish their luminosity distance-redshift 
relation, itself inferred from the flux measurement of their light curves. 

Shortly after this discovery, it was proposed by a small number of authors 
that this effect could be due to the large-scale inhomogeneities of the Universe 
[TBI [T71 ITS] . After a period of relative disaffection, this proposal experienced a 
renewed interest about five years ago. 

Now, the accelerated expansion interpretation was sufficiently misleading 
such as to induce some authors to try to derive or rule out no-go theorems, 
i.e., theorems stating that a locally defined expansion cannot be accelerating in 
inhomogeneous models satisfying the strong energy condition. But, as we have 
seen above, this is not the point. 

Other authors stressed, more accurately, that the definition of a deceleration 
parameter in an inhomogeneous model is tricky [19, 20 and has nothing to do 
with reproducing the supernova data [5T] . 

It is well-known, from the work of Mustapha, Hellaby and Ellis [32] that an 
infinite class of L-T models can fit a given set of observations isotropic around 
the observer. This has been used by Celerier [T7] to examplify her demonstration 
that models which are spherically symmetric around the observer can fit the 
supernova data and that the problem is completely degenerate. This is the 
reason why many different central observer L-T models have been proposed in 
the literature and shown to succeed rather well. 

Thus, to constrain the model further on, it is mandatory to fit it to other 
cosmological data. 



4 



density 



4 
3.5 

3 
2.5 

2 
1.5 

1 

0.5 





-40 



■■3q 2 - 



" Z0 1fl - -v<n 

X[Mpc] °10 2 - . _ .„ 10 20 30^0 

20 30 4 - 0JK) -30 -20 -10 °Y [M pc] 



30 



20 



10 



-10 



-20 



-30 



□ 



-30 -20 -10 10 20 
Y[Mpc] 



30 



1—1 3.5 

3 

2.5 
2 

1.5 
1 

0.5 




Figure 1: Double structure: the present-day density distribution, p/pb, of a 
double-structure model in back-ground units. The left panel shows the iso- 
density curves in the spatial section X-Y. The right panel shows the present day 
colour-coded density distribution in the spatial section Y-Z. White represents 
high-density and black, low-density regions 
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Figure 2: Double structure: evolution of the density profile from 100 Myr after 
the Big Bang (1) up to the present time (7). 
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Figure 3: Triple structure: the present-day colour-coded density distribution, 
p/pb, of a triple-structure model in back-ground units, showing a slice through 
the origin. 
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4.1 The direct and inverse problems 

There exist two procedures for trying to explain away dark energy with (L-T) 
models: 

1. The direct way uses a smaller number of degrees of freedom than allowed. 
Here, one first guesses the form of the parameter functions defining a class 
of models supposed to represent our Universe with no cosmological con- 
stant or dark energy, and writes the dependence of these functions in terms 
of a limited number of constant parameters. Then one fits these constant 
parameters to the observed SN la data or to the luminosity distance- 
redshift relation of the ACDM model. 

2. The inverse problem is more general. It amounts to consider the lumi- 
nosity distance Dl(z) as given by observations or by the ACDM model 
as an input and try to select a specific L-T model with zero cosmological 
constant best fitting this relation. 

Then, to avoid degeneracy, one must jump to a further step and try to 
reproduce more and possibly all the available observational data, [531 1231 

EHEZIEB]- " 

4.2 Matching observations with a single patch L— T model 
exhibiting a central void 

The L-T solutions with a central observer has been used as a first step in the 
process of reproducing cosmological inhomogencities. They model the Universe 
with the inhomogeneities smoothed out over angles around the observer whose 
location can be anywhere (this does not contradict any Copernican principle). 
This is analogous to the smoothing out over the whole space in homogeneous 
models. The use of such models must be regarded as a first approach which 
will be followed in the future by more precise ways of dealing with the observed 
inhomogeneities . 

As it has been recalled in the second section of this article, an L-T model is 
defined by two independent arbitrary functions of the radial coordinate, which 
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must be fitted to the observational data. However, in most of the proposals avail- 
able in the literature, the generality of the models have been artificially limited 
by giving the initial-data functions a handpicked algebraic form depending on 
the authors' feelings about which kind of model would best represent our Uni- 
verse. Only a few constant parameters are left arbitrary to be adapted to the 
observations. 

Another way in which the generality of the L-T models has been artificially 
limited is the assumption that the age of the Universe is everywhere the same, 
i.e. that the L-T bang-time function ts is constant. With ts being constant, 
the only single-patch L-T model that fits observations is one with a giant void 
[29l [30] . The argument brought in defense of the constant ts assumption is 
that a non-constant ts generates decreasing modes of perturbation of the metric 
[3"Tl [32] ) so any substantial inhomogeneity at the present time stemming from 
tB,r i= would imply huge perturbations of homogeneity at the last scattering. 
This, in turn, would contradict the CMB observations and the implications of 
inflationary models. 

In the recent years, we have thus seen the increase in popularity of a large, 
then huge, then giant void model, where the observer is located at or near the 
centre of a large, huge, giant L-T void of size of up to a few Gpc. The most 
achieved example of this class of models is the GBH model, specified by its 
matter content Q]\j( r ) an d its expansion rate H(r), governed by 5 free constant 
parameters and matching to an Einstein-de Sitter universe at large scales |25j . 
It is fitted to a series of observations (CMB, LSS, BAO, SN la, HST measure 
of Ho, age of the globular clusters, gas fraction in clusters, kinematic Sunyaev- 
Zel'dovich effect for 9 distant galaxy clusters) used to constrain its parameters. 
With this class of models, the conclusion is that the possibility that we leave 
close to the centre of a large (around 2.5 Gpc) void within an Einstein-de Sitter 
Universe with no dark energy is not excluded. 

However, as many other central void models, the GBH model is constructed 
with a central underdensity defined a priori by a set of constant parameters (the 
underdensity at the centre of the void, the size of the void and the transition 
width of the void profile). Hence, the outcome of the fit can be nothing but a 
central void. 

Why, then, several researchers have been led astray by the frequent claims 
of a giant void being implied by an L-T model? The reason might be twofold. 

From the results of an analysis of 44 SN la distances they performed in the 
framework of a FLRW model with £l m = 1 and fl\ = 0, Zehavi et al. [33] 
suggested a monopole in the peculiar velocity field they interpreted as marginal 
evidence for a local void of radius 70 h~ x Mpc and ^20% underdensity sur- 
rounded by a dense shell roughly coinciding with the local great walls. They 
called this putative local void the 'Hubble Bubble', and, under this name, it 
became popular in the literature. 

Strangely, a mere toy model can become nearly as strong a reference as 
actual observations. In a series of articles, Tomita [T5] used a simple toy model 
to explain in a very pedagogical way how 'dark energy' could in principle be 
mimicked by a local void. He assumed a low-density inner homogeneous region 
is connected at some redshift to an outer homogeneous region of higher density. 
Both regions decelerate, but since the inner void expands faster than the outer 
region, an apparent acceleration is experienced by the observer located inside 
this void and looking at supernovae bursting in the outer region. Such a toy 
model was only designed to stress that an actual accelerated expansion is not 
needed to reproduce the SN la data. However, many authors have taken this 
example literally and cited this model as if it were evidence for the existence of 
a local void. 
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Figure 5: The hump: the density at the present instant as a function of the 
current areal radius. 

4.3 We can do without the giant void in a central observer 
L-T model 

Since, as we show below, the arguments in favour of the limitations imposed 
above on the arbitrary functions of the radial coordinate lack sufficient strength, 
there is no need to restrict a priori their degrees of freedom. 

The arguments in favour of or against a local 'Hubble Bubble' have been 
disputed at length within the astronomical community and the latest results 
seem to point to our being located in a 'Local Sheet' (~ 7 Mpc long) which 
bounds a 'Local Void', a nearly empty region with radius at least 23 Mpc [3~4"Il3"5] . 
This seems to contradict the putative existence of a giant local void. Note, 
however (see below) that this assumed local void does not exist on our past light 
cone but in hypersurfaces t = constant and it is thus in a space-like relation to 
us. 

The arguments brought in defense of the constant ts assumption are only 
expectations that should not be treated as objective truth unless they are veri- 
fied by calculations. Such calculations have already been done, and it turns out 
that the inhomogeneities in ts needed to explain a structure of present radius 30 
Mpc is of the order of a few hundred years and that this age difference between 
the oldest and youngest region would generate CMB temperature fluctuations 
equal to AT/T = 3.44 x 10~ 6 and AT/T = -2.35 x 10~ 6 , respectively [37)] . 
This is well-hidden in the observational errors at the current level of precision. 
Therefore, there is no justification to the assumption i# = constant. 

Hence, to determine the A = L-T model best fitting the observations the 
two independent functions of the radial coordinate defining this model must be 
left totally arbitrary. We thus impose neither an hand-picked algebraic form to 
these functions, nor ts = constant [3TJ . 

Using a set of data - the angular diameter distance together with the mass 
density in redshift space - both assumed to have the same form on our past 
light cone as in the ACDM model, the two left arbitrary L-T functions are 
determined and give the mass distribution in spacetime. 

In this case, the current density profile does not exhibit a giant void but a 
giant hump and we are located in a shallow and wide basin on top of this hump 
(FIGURE 5). 



9 



Note, however, that as the giant void, the giant hump is not directly observ- 
able. It exists in the space t = now, of events simultaneous with our present 
instant in the cosmological synchronization, i.e. it is in a space-like relation 
to us. However, contrary to the giant void of the GBH type which does not 
reproduce the ACDM p(z) function and can thus be expected to be testable 
by high redshift galaxy counts, the giant hump is not observable in p(z) or in 
the number count data, since, by construction, it is designed to reproduce them 
faithfully. 

Why such a difference between the density distribution on our past light cone 
and in the t — now space? It is due to one basic feature of the L-T model (and 
in fact of all inhomogeneous models): on any initial data hypersurface, whether 
it is a light cone orat = constant space, the density and velocity distributions 
are two algebraically independent functions of the radial position. Thus the 
density on a later hypersurface may be quite different, since it depends on both 
initial functions. Whatever initial density distribution can be completely trans- 
formed by the velocity distribution. For example, as predicted by Mustapha 
and Hellaby [7] and explicitly demonstrated by Krasinski and Hellaby [9], any 
initial overdensity can evolve into a void and vice versa. In FLRW models, there 
are no physical functions of position, and all worldlines evolve together. Thus, 
while dealing with an L-T (or any inhomogeneous) model, one must forget all 
Robertson- Walker-inspired prejudices and expectations. 

It is also worth emphasizing that the existence of the giant hump is most 
probably a feature of the particular L-T models that we ended up with. This 
result must not become the starting point of a new paradigm in observational 
cosmology, aimed at detecting the hump. Before this happens, it must be de- 
cided at the theoretical level whether the hump is a necessary implication of 
L-T models properly fitted to other observations. 

4.4 An example of L—T Swiss-cheese: Marra et al.'s model 

The model is a lattice of L-T bubbles with radius >350 Mpc in an Einstein-de 
Sitter (EdS) background. Initially, the void at the center of each hole is domi- 
nated by negative curvature and a compensating overdensity matches smoothly 
the density and curvature EdS values at the border of the hole [3S a ■ 

Since the voids expand faster than the cheese, the overdense regions contract 
and become thin shells at the borders of the bubbles while underdense regions 
turn into emptier voids, eventually occupying most of the volume. 

This shows once more that the evolution of the voids bends the photon paths 
and affects more photon physics than the geometry of the inhomogeneities. 
But, here, the inhomogeneities of the model are only able to partly mimic the 
effect of dark energy. These results have induced us to consider QSS Swiss- 
cheese models which exhibit enhanced structure evolution (see the subsection 
'Structure evolution with quasi-spherical Szekeres models') and with which we 
thus hope to increase the fit to the 'dark energy' component [55] , 

5 Extracting the cosmic metric from observa- 
tions 

The inverse problem of deriving the arbitrary functions of a L-T model from 
observations is very much involved. This is the reason why most of the authors 
who have tried to deal with this issue have added some a priori constraints to 
the model [301 13T]. 
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However, Lu and Hellaby 26], McClure and Hellaby \TJ\ and Hellaby and 
Alfedeed [55] have initiated a program to extract the metric from a set of ob- 
servations. This is the full inverse problem and it is not degenerate. To date it 
has assumed the metric has the L-T form, as a relatively simple case to start 
from, though the long term intention is to remove the assumption of spherical 
symmetry. 

These authors have developed and coded an algorithm that generates the 
L-T metric functions, given observational data on the redshifts, apparent lu- 
minosities or angular diameters, number counts of galaxies, estimates for the 
absolute luminosities or true diameters and source masses, as functions of z. 
This allows both of the physical functions of an L-T model to be determined 
without any a priori assumption on their form. 

6 Conclusion 

The increasing precision of observational data implies that FLRW models must 
now be considered just a zeroth order approximation, and linear perturbation 
theory a first order approximation whose domain of validity is an early, nearly 
homogeneous Universe. 

In the nonlinear regime, which was entered since structures formed, there 
is no escape from the use of exact methods (or of averaging schemes aiming at 
investigating this issue from the standpoint of backreaction) . 

In the era of 'precision cosmology', the effect of the inhomogeneities on the 
determination of the cosmological models cannot be ignored. Inhomogeneous 
models constitute an exact perturbation of the Friedmann background and can 
reproduce it as a limit with any precision. This is the reason why they are fully 
adapted for the purpose of studying astrophysical and cosmological effects and 
for constructing precise models of universe. 

While using L-T models with a central observer to represent our 'local' 
Universe averaged over angles around us, a giant void is not mandatory to 
explain away dark energy. A giant overdensity can also do the job. However 
while neither the void nor the overdensity are directly observable, the giant void 
alone can be tested with observations of the density function on our past light 
cone. The giant hump will need more and more precise data to be constrained. 

Exact inhomogeneous solutions can be employed not only for studying the 
geometry and dynamics of the Universe, but also to investigate the formation 
and the evolution of structures. They give enhanced formation efficiency and 
might therefore help solving the problem of structure formation pertaining to 
the standard model. 

While the L-T models have been mostly used up to now for modeling the in- 
homogeneities of the Universe, the need of getting rid of spherical symmetry for 
this purpose will lead the cosmological community to consider other solutions, 
and among them QSS models, for the future developments of inhomogeneous 
cosmology. 
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